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ABSTRACT 



We present results from a sensitive search for CO 7=1 — line emission in two z > 
6.5 Lyman-a emitters (LAEs) with the Green Bank Telescope. CO 7=1 — emission 
was not detected from either object. For HCM 6 A, at 2; ~ 6.56, the lensing magnifica- 
tion factor of ~ 4.5 implies that the CO non-detection yields stringent constraints on 
the CO 7=1 — line luminosity and molecular gas mass of the LAE, L'qq < 6. 1 x 10^ x 
( Ay/300) V2 K km s-i pc^ and Mu, < 4.9 x lO'^ x (AT//300)V2 x (Xco/O.S) Mq. 
These are the strongest limits obtained so far for a z > 6 galaxy. For IOK-1, the con- 
straints are somewhat less sensitive, L'^q < 2.3 x 10^° x (A\//300)^/^ Kkm s^^ pc^ 
and < 1.9 x 10^° x {AV/3QQy/^ x (Xco/0.8) Mq. The non-detection of 
CO 7=1 — emission in HCM 6A, whose high estimated star formation rate, dust 
extinction, and lensing magnification make it one of the best high-2; LAEs for such 
a search, implies that typical z > 6 LAEs are likely to have significantly lower 
CO 7=1 — line luminosities than massive sub-mm galaxies and hyperlimiinous in- 
frared quasars at similar redshifts, due to either a significantly lower molecular gas 
content or a higher CO-to-H2 conversion factor. 

Subject headings: cosmology: observations-early universe-galaxies: evolution-galaxies: 
formation-galaxies: high-redshift-galaxies: individual (HCM 6A) 
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1. Introduction 



In recent years, multi-wavelength selection techniques have been succ essfu 



different populations of star-forming galaxies at high redshifts, z > 6 [see 



in ide ntifying 



EUisI (12008b for a 



recent review]. An important class of such galaxies consists of the Lyman-a emitters (hereafter 



LAEs), objects identified through the ir excess emiss i on in narrow 



redshifted Lyman-a wavelength (e.g. 



Hu et al, 



1998 



Rhoads et al 



jand i mages centred on the 



20001 : 



Taniguchi et al 



2005h . 



Follow-up Lyman-a s pectroscopy has yielded accurate redshifts for a significant fraction of the 



LAE population (e.g. 



Taniguchi et al. 



2005 



Kashikawa et al. 



20061) . unlike most other high-z 



star-forming galaxies (e.g. Lyman-break systems, s ub-mm ga 
photometric redshifts. In fact, an LAE at z = 6.96 (live et al. 
spectroscopic redshift of all presently-known galaxies. 

LAEs constitute a significant fraction of the star-f orming ga 



axies ), which typically only have 



2006|) has the highest confirmed 



sufficient to reionize the Universe at earlier epochs (e.g. 



Fan et al 



axy p opulation at 2; ~ 6, 



20061). The number density of 



LAEs and the typical shape of the Lyman-a emission line therefore prov ide important probes of 



physica 



Haiman 



cond i tions in the Universe ar ound the epoch of reionization (e.g. 



2002 



Kashikawa et al. 



Haiman & Spaans 



1999 : 



20061) . Equally important, the steep drop in the space density of 



quasars at 2; > 6 implies that they would be una ble to produce t he ultraviolet (UV) background 



Fan et al. 



20011). Star-fo rming galaxies like the 



Yan & Windhorst 



2004) 



radiation required to reionize the Universe (e.g 

LAEs are thus most likely to have been responsible for reionization (e.g 
Understanding the factors that influence the star-formation activity in these systems (e.g. the 
molecular gas content, the star formation efficiency, etc) is hence of much importance. 

Recent studies of individual LAEs at 4.2 < z < 6.6 have shown that these galaxies have a 



3 - 1000 Myr and m etallicities from 0.005 - 1 Zq (e.g. 



Charv et al. 


2005; 


Lai et al. 


2007 



Finkelstein et al 



20091). Evidence has also been found for significant amounts of dust in LAEs, 
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with UV extinctions of A 



Chary et al. 



z > 6 (e.g. 



UV 



0.5 - 5 mag. (e.g. ICharv et al.ll2005l : Finkelstein et all 120091) . 



(l2005h suggest that the observed decrease in the the cosmic star- formation rate at 



Stanwav et al. 



2004|) might be explained by the presence of dust in the z > 6 LAE 



population. Typical star-formation rates (SF Rs) in LAEs, as inferre d from their UV continua, are 



low, a few tens of solar masses per year (e.g. 



Taniguchi et al. 



2005h : these are 2 — 5 times larger 



than the SFRs derived from the resonantly- scattered Lyman-a line. While this difference could 
arise from dust obscuration effects, Hi absorption by the damping wing of the Gunn-Petersen 
trough is likely to also contribute towards reducing the strength of the Lyman-o; line for LAEs at 



z > 6 .5 (iHaiman 



LAE (Chary et al. 



20021). We note that the tentative detection of the Ha line in a single 2; ~ 6.56 
2005|) yielded a significantly higher SFR estimate (140 Mq yr^^) than that 



obtained even from the rest-frame UV continuum (~ 9 Mq yr ^); this emphasizes the possibility 
that the SFRs in other LAEs might have been under-estimated due to dust extinction. 

The detectability of high-2; galaxies like the LAEs relies on their undergoing an elevated 
level of star-formation activity, which naturally requires fuel in the form of molecular gas. 
Such gas is most effectively stu died through observations of redshifted CO emission lines (e.g. 



Solomon & Vanden Bout 



20051) . The luminosity in the low-7 CO lines can be used to estimate 
the total molecular gas mass fueling the star-formation activity, while the CO line widths provide 
a measure of the dynamical mass of the galaxy. Studies of molecular gas at high redshifts, z > A, 
have so far f ocused on the most massive, far-infrared-lu minous systems, the sub-mm galaxies or 



quasars (e.g. 



Schinnerer et al. 



2008 



Walter et al 



20031), and no information is available in the 



literature on the molecular gas content of "normal" star-forming galaxies, such as the LAEs. 

In the present work, we address this outstanding problem by conducting a sensitive search 
for low-excitation CO line emission in two z > 6.5 LAEs with the Green Bank Telescope (GBT), 
which allows us to place strong constraints on their total molecular gas masses. Throughout 



this Letter, we adopt a cosmological model with (i^A, ^m, h) = (0.73, 0.27, 0.71) (|Spergel et al. 
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Fig. 1.— GBT spectra in the redshifted CO 7=1 - line from HCM 6A (left panels, [A] and 
[C]) and IOK-1 (right panels, [B] and [D]). The spectra in the top and bottom panels have been 
smoothed to velocity resolutions of 50 km s^^ and 150 km s^^, respectively. The velocity scale of 
each spectrum is relative to the Lya redshift, while the error bar in each panel indicates the redshift 
uncertainty (conservatively, Az = 0.01). 



2007 ). 



2. The targets 



Our two target LAEs were selected to have z > 6.5, to ensure that their redshifted CO 7=1 — 
lines could be observed with the sensitive GBT Ku-band receiver. The first system, HCM 6A at 
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z ~ 6.56 (IHu et al.ll2002l) . is perhaps the best candidate LAE for a searc h for molecular em ission 



as it is strongly lensed by a foreground cluster (Abell 370, at 2; ~ 0.375; 



Kneib et al 



1993), with a 



magnification factor of ~ 4.5; this significantly im proves the sensitivi ty to CO emission. It also 



has the highest estimated SFR of all known LAEs; 



Chary et al. 



( 2005 ) found excess broad-band 



emission in a Spitzer IRAC 4.5/im image, compared to the continuum at other wavelengths, 
and argue that this is due to strong Ha emission, with an implied SFR of 140 Mq yr^^. The 
SFR inferred from the UV continu um (uncorrected for dust effects) is significantly lower than 



this. ~ 9 M.^ vr~^ (IHu et al 



Taniguchi et al. 



20021), similar to values obtained in other LAEs at z ~ 6.5 (e.g. 



20051 ). If the SFR derived from the Spitzer image is correct, the discrepancy 
between the two SFR values implies a high dust extinction. Any ~ 2.6 ma g, consistent with the 



valu e independently derived from a fit to the broad-band photometric data (|Chary et al. 



also 



Schaerer & Pelld 



2005 



; see 



infrared galaxies (U 



Downes & Solomon 



20051). Such a high SFR would be similar to that of nearby ultra-luminous 
IRG s), which hav e typical molecular gas masses of 10^^^° Mq (e.g. 



19981) . Conversely, 



Boone et al 



(120071) argue that the limits on the 850 fim 



and 1.2 mm continuum flux densities of HCM 6A imply upper limits to the SFR in the range 
10 — 90 Mq yr^\ for assumed dust temperatur es in the range 18 — 54 K; these SFR estimates 



are somewhat low er than the values obtained by 



Chary et al 



Chary et al. 



(l2005b and 



Schaerer &Pell6l ( 2005 ) 



(|2005[) argue that the properties of HCM 6 A are similar to those of other z >Q LAEs; 
the molecular gas properties derived for this object are thus likely to be representative of the entire 
population. 



Our second target, IOK-1, was discovered in a Subaru survey for 2; ~ 7 LAEs (live et al. 



2006h and has the highest spectroscopically-confirmed redshift of all known galaxies. Jye et al. 



(120061) obtain a star-fo rmation rate of (10 ± 2) Mq yr ^ from the Lyman-a line luminosity. 



Taniguchi et al 



(|2005l) find that the SFR estimated from the UV continua of z ^ 6.5 LAEs is 
typically ~ 5 times larger than that inferred from the Lyman-a line, implying that the SFR of 
IOK-1 is likely to be > 50 Mq yr^^, among the highest of the z > 6 LAE population. This 
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system has not so far been observed at sub-mm or infrared wavelengths and no information is thus 
available on its dust properties. 



3. Observations and data analysis 

The search for CO 7=1 — emission {I'rcst = 115.2712 GHz) from the two LAEs was carried 
out with the GBT 12.0 — 15.4 GHz Ku-band receiver between August and September 2008, 
during excellent summer observing conditions. The Auto-Correlation Spectrometer was used as 
the backend, with two circular polarizations and a bandwidth of 200 MHz, sub-divided into 8192 
channels and centred on the redshifted CO 7=1 — line frequencies (14.48 GHz and 15.25 GHz). 
This yielded a total velocity coverage of ~ 4000 km s^^ and an initial velocity resolution of 
~ 1 km s~^, after Hanning-smoothing and re-sampling. Dual-beam nodding was used to calibrate 
the system bandpass, with a nodding timescale of 2 minutes. The telescope pointing was corrected 
every 2 hours by observations of nearby bright calibrators. System temperatures were measured 
online by firing a noise diode, and were typically found to lie in the range 22 — 30 K. The 
standard flux calibrators, 3C48, 3C147 and/or 3C286 were used to verify the absolute flux scales 
at the observing frequencies; we estimate that the uncertainty in the flux calibration is less than 
10%. The total observing times were ~ 25 hours for HCM 6A and ~ 23 hours for IOK-1, 
including overheads associated with dual-beam nodding, pointing and flux calibration. Note that 
the fuU-width-at-half-maximum of the GBT Ku-band beam is ~ 50" (corresponding to a spatial 
scale of ~ 280 kpc at z = 6.56), implying that all CO 7=1 — emission from the LAEs (and any 
nearby galaxies at similar redshifts) would lie within the telescope beam. 

The data were analyzed using the GBTIDlJU data analysis package, following standard 
procedures. A second-order baseline was fit to each scan during the initial calibration, after which 



' http ://gbtidl. sourceforge.net 
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Table 1: Results from the GBT observations of HCM 6A and IOK-1. 



Source ziya t'coi-o On-source time RMS* L'qq (Sa) (Sa) 

GHz Hrs. /iJy x(Ar/300)i/2 x (AV/300)i/2(Xco/0.8) M( 

K km s^^ pc^ 

HCM 6A 6.56 15.25 12 22.0 < 6.1 x 10^ t < 4.9 x 10^ t 

IOK-1 6.96 14.48 11 17.3 < 2.3 x 10^° < 1.9 x 10^° 

*The quoted values of RMS noise, L'qq and are at a velocity resolution of 300 km s"^. 

^For HCM 6A, the luminosity and mass limits have been corrected for a magnification factor of 

4.5. 

the data were examined on a scan-by-scan basis for residual baseline structure; this was done 
independently for both polarizations. Conservatively, all scans showing baseline structure on 
scales of a few tens of MHz (that might mimic a putative spectral line) were edited out and not 
included in the final analysis; this resulted in the excision of ~ 25% of data for each target. The 
final on-source integration times were ~ 12 hours for HCM 6 A and ~ 11 hours for IOK-1. For 
each source, data from the two polarizations were combined together with appropriate weights, 
based on their root-mean-square (RMS) noise values to produce the final spectra. The spectra 
were then smoothed to coarser resolutions (50, 150 and 300 km s^^) to search for CO 7=1 — 
line emission. 

4. Results 

The final spectra obtained towards HCM 6A and IOK-1 are shown in the four panels of 
Fig.fll at velocity resolutions of ~ 50 km s^^ and ~ 150 km s^^. The final RMS noise values per 
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50 km channel are 52 yuJy (HCM 6A) and 51 /iJy (IOK-1), and per 150 km s"^ channel are 
31/iJy (HCM 6 A) and 30/iJy (IOK-1). No evidence for CO 7=1 — line emission was seen in the 
spectra at these (or other) velocity resolutions. 

The non-detections of CO 7=1 — line emission place strong constraints on the CO line 
luminosities, which can, i n turn, be used to derive li r nits on the total mass of cold molecular gas in 



the two LAEs. Following 
can be written as 



Solomon & Vanden BoutI dlOOSb . the CO 7=1 - line luminosity L'co 



L' 



3.25 X 10^ Sco u-? Dl (1 + zV^ 



^CO - ^-^o Lyj uco '^obs -^L -T ^; , (1) 

where z/obs is the observing frequency (in GHz), is the luminosity distance (in Mpc), and L'qq 
is in K km s^^ pc^. For our non-detections, Sco AV" = 3crAV' x (in Jy km s~^) gives the 
3cr upper limit on the integrated flux density in the CO 7=1 — line, where ctav is the RMS 
noise at the velocity resolution AV. We will assume a line width of = 300 km s^\ similar 



CariUi & Wang 



20061) . We then obtain 3a 



to the median observed width in high- 2; quasars (e.g. 
CO 7=1 - line luminosity limits of L'^q < 6.1 x 10^ x (AV/300)i/2 ^ km s^^ pc^ (HCM 6A) 
and L'co < 2.3 x 10^° x (Ar/300)i/2 ^ km s^^ pc ^ (IOK-1), after co rrecting the HCM 6A 



luminosity for the lensing magnification factor of 4.5 (IKneib et al. 



1993 ) 



To convert the CO 7=1 — line luminosity to an estimate of the total molecular 



Solomon & Vanden Bout 



gas 



20051) . For 



mass requires the CO-to-H2 conversion factor Xco (e.g. 
virialized molecular clouds in a quiescent galax y like the Milky Way, the conversion factor is 
typically Xqo ~ 4.6 Mq (K km s^^ pc^)^^ (e.g. ISolomon & Barrettill99lh : conversely, luminous 
infra red galaxies have far lower CO-to-H2 conversion factors [Xqo ~ 0.8 Mq ( K km s~^ pc ^)~^; 



Downes & Solomon 



Weifi et al, 



1998n . For the nearby dwarf starburst galaxy, M82, 
(1200 It) find that the lowest values of the conversion factor are measured towards the central 
star-forming regions where the UV radiation field is most intense. Given that the two LAEs 
observed here appear to be undergoing elevated levels of star-formation, we adopt the conversion 
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factor Xco = 0.8 Mq (K km s pc ) . N ote that this conve rsion factor is also used for CO 



line studies of high-z sub-mm galaxies (e.g. 



GreveetaL 



20051), allowing a direct comparison 



between the inferred molecular masses of sub-mm galaxies and LAEs. Using this conversion 
factor then yields Mr^ < 4.9 x 10^ x (Ar/SOO)^/^ ^ (Xco/0.8) Mq for HCM 6A, and 

< 1.9 X 10^° X (Ay/300)^/2 X (Xco/0.8) Mq forIOK-1. We emphasize that the upper 
limits to the molecular gas mass would be higher than the above values if the Milky Way 
conversion factor were applicable to these LAEs. The results of our observations are summarized 
in Table [B 

Finally, the CO 7=1 — line luminosity and the FIR luminosity are correlated in nearby 



([Riechers et al 



starburst and spiral ga laxies (IGao & Solomonll2004l) . with logLpTR = (1.26±0.08) xlogLcQ — 0.81 



20061). Our limits on the CO 7=1 — line luminosity can be combined with 
this relation to infer the FIR luminosity of the two LAEs and thence, their SFRs. We obtain 

Lfir < 3.3 X 10^1 Lq (HCM 6A) and Lpm < 1.8 x lO^^ (IOK-1), and SFRs of < 66 Mq yr~^ 
(HCM 6A) and < 360Mq yr'^ (IOK-1), on using the relation SFR = 2 x 10-^^{Lfi^/Lq)Mq yr"^ 



(Kennicutt 



19981). The limit to the SFR in HCM 6A is similar to that obtained from the 1.2 mm 



and 850 /xm continuum imaging (SFR < 10 — 90 Mq yr ^; 



Boone et al. 



2007 ). While our SFR 



limit in HCM 6A is signi ficantly lower than the estimate of 140 Mq yr from the tentative 



detection of Ha emission (jChary et al. 



I2005 I). it should be pointed out that our result depends on 
the assumption that the local correlation between CO 7=1 — line luminosity and FIR luminosity 
is applicable in 2; > 6 LAEs. Note, however, that some classes of high-redshif t galaxies [e.g- the 



"B^K" galaxies, selected as outliers in plots of (B — z) versus (z — K) colors; 



Daddi et al. 



200411 



have larger CO line luminosities than predicte d by the FIR-CO r elation, with large molecular 



gas masses but low star formation efficiencies (IDaddi et al. 



20081). If the high-redshift LAEs are 



similar to the BzK galaxies, our limits to the CO 7=1 — line luminosities would imply even 
lower FIR luminosities, and star-formation rates, than those listed above. 
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5. Discussion 



These are the first constraints on the CO 7=1 — line luminosity and the molecular gas mass 
of 2 > 6 LAEs, providing a new window into physical conditions in star-forming galaxies at the 
highest redshifts. The high inferred SFR and dust extinction in HCM 6A, as well as its large 
magnification factor, imply that it is one of the best candidates for a detection of CO emission 
in a high-z LAE. The limit on its CO 7=1 — line luminosity obtained here is the deepest ever 
obtained for a 2; > 6 galaxy. For the assumed CO-to-H2 conversion factor, the m olecular gas ma ss 



limit is within a factor of ~ 2 of the molecular gas mass of the Milky Way (e.g. 



Combes 



1991 ) 



While the limit on the CO 7=1 — line luminosity of IOK-1 is not as strong, i t is lower than the 



median luminosity in high-2; sub-mm galaxies (T^co ~ 3.8 x 10^° K km s ^ pc^; 



Greve et al. 



2005 ) 



and similar to the luminosities of lower-redshift BzK galaxies (which have similar SFRs; e.g. 



Daddi et al. 



20081). Our results for both LAEs rule out the presence of extreme CO 7= 1 — line 



luminosities, such as tho s e seen in sub-mm g alaxies or hyperluminous IR quasars (e.g. 



Greve et al. 



2005 



Walter et al. 



2003 



Carilli et al 



20071). LAEs thus appear to either contain significantly 



lower quantities of cold molecular gas or have significantly higher CO-to-H2 conversion factors 
than sub-mm galaxies or FIR-bright quasar hosts. The latter possib ility cannot be ruled out as the 



19881), and 



CO-to-H2 conversion factor is likely to depend on metallicity (e.g. Maloney & Black 
quasar host galaxies and sub-mm galaxies appear to be dusty, metal-rich systems. 

Finally, it is clear that molecular gas must be present in LAEs to fuel the observed 
star-formation activity and dust reddening. In cases where the CO line emission is optically thick 
and thermalized, the flux density in the CO lines scales oc z/^ oc Ju, where J\j is the rotational 
quantum number of the upper level. This suggests that, despite the sensitive limits obtained 
here, the Ju > 5 CO lines may provide a more effective avenue to probe the molecular gas 
content of high- 5; LAEs, with planned facilities like the Expanded Very Large Array (EVLA) 
and the Atacama Large Millimetre Array (ALMA). For example, for optically-thick, thermalized 
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emission, ALMA would be able to detect the CO 7=6-5 line from a z = 6.5 star- forming galaxy 
with =3x10^ Mq in 3 hours of on-source integration time. Unfortunately, the high kinetic 
temperatures and densities required to raise the CO molecules to the high- J excitation states are 
unlikely t o be p resent in "normal" star-forming galaxies like the LAEs. For example, Fig. 9 of 



Aoetal 



(|2008h shows that the CO line intensities are sub-thermal at Ju > 5 in almost all galaxies 
(including ULIRGs and sub-mm galaxies) with observations of these lines. The sole exception 
is the lensed quasar APM08279-I-5255, where a c ombination of AG N heating and very high gas 



densities appears to yield the high CO excitation (IWeiB et al 



20071) . This implies that it is likely 



to be difficult to detect the high- J CO lines from z > 6 LAEs even with ALMA. The 158;um 
fine- structure transition of ionized carbon may hence prove the best candidate for mapping the 
arge-scale structure of high- 2; star-forming galaxies and determining their dynamical masses 



( Walter & Carilli 



20081). A search for this transition in HCM 6 A is currently in progress [Kanekar 



et al. (in prep.)]. 



6. Summary 

We have carried out a deep GBT Ku-band search for CO 7=1 — line emission from two 
high-redshift Lyman-a emitters, HCM 6A at z ~ 6.56 and IOK-1 at 2; ~ 6.96. Our non-detection 
of CO emission from the lensed LAE, HCM 6A, implies strong constraints on its CO 7=1 — 
line luminosity and molecular gas mass, L'^q < 6.1 x 10^ x (Ay/300)^/^ K km s^^ pc^ 
and < 4.9 x 10^ x (AF/300)i/2 ^ (Xco/0.8) Mq, the strongest obtained to date 
for a 2; > 6 galaxy. In the case of IOK-1, the absence of a lensing magnification factor 
implies that the limits on the line luminosity and molecular gas mass are somewhat less 
sensitive than for HCM 6A; we obtain L'^o < 2.3 x 10^ x (Ay/300)i/2 ^ km s^^ pc^ and 
< 1.9 x 10^° X (AF/300)i/2 X (Xco/0.8) Mq. If HCM 6A is a typical LAE at these 
redshifts, our results imply that LAEs are unlikely to show high CO 7=1 — line luminosities 
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(such as those found in quasar host galaxies at similar redshifts), due to either a lower molecular 
gas content or a higher value of the CO-to-H2 conversion factor. It hence appears that observations 
of redshifted CO line emission are unlikely to be an effective means of studying gas dynamics in 
the less-luminous galaxies responsible for reionizing the Universe, even with upcoming arrays 
such as ALMA and the EVLA. 
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